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Summary 


/frhis  paper  develops  a  theoretical  analysis  to  determine  the  force  level  attained 

/jll 

in  a  cylindrical  metal  tube  when  the  _Ube  is  fragmenting  or  splitting  and  rolling 
on  a  rigid  die*  The  tube  is  treated  as  a  thin,  rigid— perfectly  plastic  cylindrical 
shell  which  fractures  when  a  critical  level  of  maximum  strain  is  attained.  The  die 


has  a  semi* 


•toroidal  shape  in  order  to  expand  the  tube  radially  and  to  bend  it  meri- 


dianally,  thus  leading  to  fragmentation,  i,e«,  fractures  running  in  two  perpendicu- 

/)/i 

lar  directions^/lt  is  shown  that  in  the  fragmentation  regime  the  meridianally  run¬ 
ning  cracks  advance  rapidly  at  first  and  then  may  become  arrested.  It  is^  also  shown 
that  localized  high  bending^ strains  develop  and  lead  to  transverse  fracture/}  The 
energy  absorbed  during  the  process  is  evaluated  through  a  numerical  procedure  which 
considers  successive  configurations  of  the  deforming  tube.  The  effect  of  contact 


friction  on  the  energy  absorbed  is  considered  and  found  to  be  significant. 


J4lA 
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Nomenclature 


The  tube  mean  radius  and  thickness  are  A  and  T,  respectively.  All  dimensionless 

length  parameters  are  normalized  with  respect  to  the  shell  parameter 

a  A/^'/aT,  dimensionless  tube  mean  radius 

b  see  eq.  (25) 

D  dimensionless  displacement  of  the  tube 

f  see  eq.  (12) 

F  die  contact  force  per  unit  circumference 

g  ^  -  coordinate  of  the  plastic  hinge 

h  dimensionless  tube  half  thickness 

H(s)  Heaviside  unit  step  function  of  s 

Z  dimensionless  petal  length, 

Z  lower  limit  of  integrals,  see  eq.  (l5) 

m  M/cJqT  ,  dimensionless  bending  moment 

M  bending  moment  per  unit  length 

n^  N^/aoT,  dimensionless  meridianal  force 

meridianal  force  per  unit  length 
n  N  /anT,  dimensionless  circumferential  force 

circumferential  force  per  unit  length 
p  dimensionless  length  of  plastic  region, 

p  contact  pressure 

q  ^  dimensionless  shearing  force 

Q  shearing  force  per  unit  length 

R  dimensionless  die  radius 

dimensionless  modified  die  radius 
t  2h  «  T/n/at',  dimensionless  tube  thickness 

u  dimensionless  radial  displacement  of  middle  surface 

u'  dimensionless  radial  displacement  of  inner  surface 

V  dimensionless  radial  velocity 

V  radial  velocity 

w  ^  dimensionless  tube  axial  velocity 

d^ 

X  dimensional  position  coordinate 

y  dimensionless  length  projection  on  tube  axis 

z  dimensionless  projection  of  petal  on  tube  axis 

6(s)  Dirac  delta  function  of  s 


iv 


t 

& 

•n 


V 

X 

o 

K 

X.,X^  A 

V 

C 


a 


o 


CO 


average  circumferential  strain 
circumferential  strain  rate 

x/'vSt*,  dimensionless  position  coordinate  measured  from  the  plastic 
hinge  in  the  shell 

dimensionless  curvature 
dimensionless  curvature  rate 
dimensional  curvature  rate 
arbitrary  positive  multipliers 
friction  coefficient 
angle j  see  eqo  (26) 

dimensionless  position  coordinate  measured  from  the  petal  tip 

compressive  stress  in  the  undeformed  portion  of  the  tube*  Prime  indicates 
that  the  effect  of  axial  force  on  the  yield  surface  has  been  considered. 

material  flow  stress  in  tension  or  compression 

dimensionless  time  parameter 

slope  angle 

function  defining  flow  surface 
die  angle 

dimensionless  angular  velocity 


Subscripts 

c  refer  to  location  at  crack  tip  or  petal  root 

h  refer  to  location  at  plastic  hinge 

t  refer  to  location  at  petal  tip 

o  initial  condition,  except  for  o 


V 


Introduction 


A  method  proposed  for  the  absorption  of  kinetic  energy  is  the  process  of  frag¬ 
mentation  of  metal  tubing,  A  metal  circular  tube  is  forced  over  a  die  causing 
circumferential  expansion  and  meridianal  bending  of  the  tube  as  shown  in  figures  1  (^) 
and  1(b),  Due  to  the  limited  ductility  of  the  material, cracks  develop  along  meri¬ 
dianal  lines  so  that  a  number  of  petals  are  formed  at  the  end  of  the  tube.  In  the 
fragmenting  regime  these  petals  grow  and  then  break  off  due  to  the  large  bending 
strains  developed.  The  petal  growth  and  breaking  process  then  repeats.  After  a  few 
petals  have  broken  off,  the  process  becomes  irregular  so  that  at  a  given  time  the 
petals  have  attained  different  stages  of  development.  Energy  is  absorbed  essen¬ 
tially  by  plastic  flow  of  the  tubing  material  and  by  contact  friction.  The  energy 
absorbed  by  the  fracture  process  is  small  in  comparison  with  that  absorbed  by 
plastic  flow  and  friction. 

The  purpose  of  this  paper  is  to  develop  a  theoretical  analysis  to  determine  the 
force  level  in  the  tube  and  whether  or  not  the  petals  will  break  off,  A  quasi¬ 
static  process  is  assumed, i,e,,  inertia  terms  are  neglected  in  the  equations  of 
motion.  This  means  that  the  analysis  is  appropriate  only  when  the  maximum  kinetic 
energy  associated  with  the  deformation  of  the  tube  is  small  in  comparison  with  the 
total  energy  absorbed  and  the  velocities  are  small  in  comparison  with  the  elastic 
wave  propagation  speed. 

McGehee  [  l]^  has  conducted  experiments  on  the  fragmentation  of  202U-T3 
aluminum  alloy  tubing. 


lumbers  in  brackets  refer  to  references  listed  at  the  end  of  the  paper 


II  Analysis  of  Petal  Growth 


Formulation  of  the  Problem 

The  theoretical  analysis  of  the  growth  of  petals  or  cracks  is  based  on  seven 
assumptions o  These  assumptions  follow* 

a  At  a  given  instant  all  petals  have  the  same  configuration.  This  assumption  is  appro- 
priate  if  one  wishes  to  determine  the  force  in  the  tube  due  to  deformation  of  one  petal 
or  the  average  force  in  the  tube.  The  theoretical  peak  tube  stresses  found  from  the 
following  theoretical  analysis  occur  only  before  the  irregular  process  develops, 
b  The  propagation  of  meridianal  cracks  is  governed  by  a  condition  of  critical  radial 
displacement  at  the  crack  tips^  hence,  the  average  circumferential  strain  along  a  circle 
passing  through  the  meridianal  crack  tips  is  an  empirical  constant. 

c  The  loading  and  deformation  of  the  unsplit  portion  of  the  tube  is  axisymmetric.  This 
assumption  follows  from  a  provided  a  sufficiently  large  number,  say  six  or  more,  petals 
form. 

d  The  assumptions  of  the  theory  of  thin  shells  apply  to  the  unsplit  portion  of  the  tube. 
This  is  probably  the  least  accurate  assumption  made,  since  the  tube  thickness  to  radius 
ratios  employed  in  [  l]  are  not  sufficiently  small.  However,  without  this  assumption 
an  analysis  would  be  virtually  impossible. 

e  The  tube  is  homogeneous  and  is  made  of  a  rigid-perfectly  plastic  material  which  obeys  the 
Tresca  or  maximum  shear  stress  yield  condition  and  associated  flow  rule.  This  assun^jtion 
is  reasonable  provided  the  plastic  strains  are  large  in  comparison  with  the  elastic  strain 
at  yielding  and  the  material  has  negligible  work  hardening  and  anisotropy, 
f  Geometry  changes  are  small  enough  so  that  equations  of  equilibrium  and  compatibility  of 
the  unsplit  portion  of  the  tube  are  those  of  a  cylindrical  shell, 
g  Coulomb  friction  occurs  at  all  points  of  contact. 

In  addition,  the  analysis  requires  two  a  priori  conditions.  These  follow, 
h  Contact  between  the  tube  material  and  the  die  occurs  only  at  the  tips  of  the  formed 
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petals o 

1  The  current  meridianal  curvature  of  the  petals  was  produced  when  that  point 
currently  in  the  petal  was  part  of  the  plastically  deforming  shell.  The  petals  are 
assumed  to  be  rigid^i.e.,  have  unloaded  from  a  plastic  state. 

Figure  2(a)  shows  the  cross  section  of  the  tube  CHI  with  a  petal  TC  at  its  end. 
Contact  with  the  die  occurs  at  point  T’ .  The  region  CHI  behaves  as  a  semi  infinite 
plastic  cylindrical  shell  subjected  to  axisymmetric  transverse  shear  and  bending  moment 
at  Go  For  a  development  of  the  plastic  analysis  of  such  a  shell  refer  to  the  Appendix^ 
Figure  2(b)  shows  the  free  body  diagram  of  the  petal  and  tube. 

The  dimensionless  lengths  appearing  in  the  analysis  are  normalized  with  respect 
to  the  shell  parameter  nTaT.  The  dimensionless  stress  resultants  and  m 

are  defined  in  the  Appendix  and  in  the  Nomenclature.  5  is  a  dimensionless  position 
coordinate  measured  from  the  petal  tip.  The  lengths  ^  “  TC/'^/aT,  p  “  CH/JaT,  and 

g  =  +  P,  (1) 

shown  in  figure  2(a)  all  change  with  time.  The  dimensionless  radial  displacement  u, 
the  slope  angle  (p,  and  the  dimensionless  meridianal  curvature  x  are  functions  of 
the  position  coordinate  ^  and  time. *  The  subscripts  t,  c,  and  h  refer  to  loca¬ 
tions  T,  C,  and  H,  respectively. 

Let  u  =  u(?^t),  where  t  increases  monotonically  with  the  time.  Since  the 
material  properties  are  time  independent^  any  variable  which  increases  monotoni¬ 
cally  with  time  may  replace  't.  From  assumption  b 

u(^,'r)  a  constant  (2) 

where  is  the  average  circumferential  strain  at  ^  ^ 

From  differentiating  (2) 


9^ 


0 


or 
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The  slope  angle  at  a  given  material  point  in  the  plastic  region  is  the  sum 
of  the  initial  slope  angle  and  the  additional  slope  developed  by  plastic 
flow  during  the  crack  growth  process* 

Thus 


<P  (Cj'T)  =  <Pq(5)  +  Jf^  d  P  (U) 

As  is  developed  in  the  Appendix,  the  dimensionless  velocity  v,  angular  velocity 
(X),  and  curvature  rate  x  depend  on  the  dimensionless  position  coordinate  t]  measured 
from  the  plastic  hinge  and  an  arbitrary  multiplier.  At  a  given  instant  is  pro¬ 
portional  to  u)(r|)  where  “n  =  g  -  ?•  Also,  is  in  the  same  proportion  to  v(r))* 

dr 

Then  it  follows  that 


(5) 


Now  let  T  ~  .  The  solutions  show  that  S  increases  monotonically  with  timej 

hence,  this  step  is  admissible.  Substitute  (5)  and  (3)  into  (U),  and  let  ^ 

to  obtain 

=‘Po(^)  *  f  w[  .g(s)-  i]  <p(s)ds  (6) 

vlp(s)]  ° 


Prom  the  relationship 


(7) 


one  finds  an  integral  expression  for 

X^(^)  =  x^(^)  +  j  dfifl^<Pc(s)ds  (8) 

^  o 

Since  the  petals  remain  rigid,  the  curvature  of  the  petals  is 

x(C)  =  Xj.(C)  (5') 

Into  (6)  and  (8)  are  substituted  expressions  for  the  spacial  distributions  of 
the  dimensionless  radial  velocity,  angular  velocity,  and  curvature  rate  in  a  rigid- 
plastic  shell  subjected  to  axisymmetric  bending  moment  and  shear  force  at  its  end.  In  the 
equations  which  follow  the  approximate  expressions  forv,  co,  and  x  (equations  (83), 
(8H),  and  (85),  respectively,  of  the  Appendix)  are  used  in  place  of  the  exact  ex¬ 
pressions  derived  in  the  Appendix.  This  produces  some  sinplification  at  the  ex¬ 
pense  of  only  slight  errors.  Since  the  variable  p  is  always  positive,  one  may  omit 
the  Heaviside  function  as  a  factor  in  v[p(s)].  One  then  obtains 


-  fVIo  f 


f(s).H[g(s)-^]ds 


•f^/o  f(s)*g(s)*H[g(s)  -«]ds 

A 

X(.(^)  =x^(-«)  +  J  f(s)*6[g(s)  -  ^]ds 

f(s)*H[g(s)  -  i]ds 
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where 


f(s)  = 


«Pc(s)‘H(s) 
p(s)  +  ^-y/IlpCs)] 


The  definite  integrals  having  the  Heaviside  step  function  and  the  Dirac  delta 
function  as  factors  of  the  integrands  are  evaluated  by  the  following  formulas: 


a  <  b 


f(x)*H[g(x)  -  c]±x: 


r 

f(x)dx,  c<  g(a) 

? 

I  f (x)dx,  g(a)  <  c  <  g(b) 
g"^(c) 

0  ,  c  >  g(b) 


0  ,  c  >  g(a) 


f(x)»  6[g(x)  -  c]dx 


f[g~^(c)]^|^  >g(a)  <  c<  g(b) 

^  g=c 


0  ,  c  >  g(b) 


where  g”^  is  the  inverse  of  the  function  g,  i.e.  if  y  =  g(x),  then  x  =  g"^(y). 


Equations  (lO)  and  (ll)  become 


<p^(^)  =<Po(<)  +  (1  -1^3  f (s)ds 


f(s)*g(s)ds 


xU) 

C 


xJS)  *  f(S^) 


(2h) 


where,  with  the  initial  value  of  p,  the  lower  limit  i*  satisfies 

S*  ‘  0  ^  <  Po 

I*  =  g~^{z)  ^  >  Po 


(15) 


The  ^  diagram  of  figure  3  aids  the  understanding  of  equations  (13) >  (lU)j  and 

(15).  The  line  OT  represents  the  5 -coordinate  of  the  petal  tip,  OA  that  of  the  crack 
tip,  and  the  plastic  hinge  in  the  shell  as  I  increases  from  zero.  The 

region  to  the  left  of  OA  is  the  petal.  That  between  OA  and  is  the  plastic 

region  in  the  shell,  and  that  to  the  right  of  H°H'“h®HA  is  rigid.  Lengths  shown  in 
the  diagram  are 

OT  =  TO  =  i  ,  TH  «  g{i),  CH  =  g(^)  -  ^  =  p(^) 


OT®  =  s,  T®H®  =  g(s),  BH®  =  g(s)  -  i  >  0 

OT*  ”  I*,  T*H*  •=  g(^*)  =  TO  = 


0H°  =  OT®  =  Po  =  go 


A  vertical  line  in  the  diagram  represents  states  through  which  a  material  point 
must  pass«  If  the  crack  tips  have  not  reached  the  point  then  the  material 
point  currently  at  the  crack  tip  was  plastic  at  the  time  when  the  crack  began  to  grow. 
In  this  case  plastic  deformation  at  the  material  point  occurs  during  the  time  interval 
from  0  to  ^  •  Thus, the  lower  limit  of  the  integral  is  =  0,  After  the  crack  tips 
have  passed  point  to  point  C,  the  material  point  currently  at  C  was  initially  rigid. 
It  underwent  no  plastic  deformation  until  the  plastic  hinge  reached  that  point.  Thus, 


7 


the  lower  limit  of  the  integrals  must  be  the  crack  length  at  the  instant  the  point 
currently  at  C  became  plastic.  In  the  diagram  this  length  is  =>  OT*,  found  by 
dropping  a  vertical  line  from  C  to  and  then  drawing  the  horizontal  line 
Since  the  equation  of  the  line  is  €=  gC-^)  or  ^  =  g"^(C)j,  then 

The  shape  of  curve  shown  is  typical  for  many  confuted  cases. 

As  it  is  shown  in  the  Appendix,  the  dimensionless  length  of  the  plastic  region 
depends  on  the  ratio  of  the  applied  bending  moment  to  applied  shearing  focus.  This 
ratio  depends  on  the  instanteous  geometry  and  rate  of  change  of  geometry  of  the  petal. 
The  following  shows  how  to  find  the  instantaneous  configuration. 

From  the  condition  that  the  petal  is  rigid  one  finds  the  slope  angle  <p(S,^)  at 
any  point  in  the  petal  to  be 


(p(^/)  =  (p  +  r  X  (?)d5  (16) 

C  J  ^ 

K 

The  tip  slope  angle  tip  radial  displacement  U|.(^)  and  the  projection  z  ) 

of  the  petal  center  line  (TC  of  figure  2(a))  on  the  tube  axis  are 


(17) 

(18) 

(19) 


Since  a  numerical  technique  will  eventually  be  used  to  obtain  solutions,  it  is  more 

convenient  to  put  equations  (17),  (18),  and  (19)  in  derivative  forms  rather  than 

integral  forms.  Taking  the  derivative  of  <P^  with  respect  to  ^  and  substituting 

(13)  and  (lli)  for  cp  aiid  x  results  in 
c  c 
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(20) 


=  [1  +  |^i3  p(^)3  f(^) 

By  taking  the  derivative  of  and  zt  with  respect  to  I  and  noting  that  for  a 
rigid  petal 

d  d  <p. 

=  -dT- 


one  finds 


du^  d(^ 

—  =  sin<p^  H  dT 


dzt  <i  <P^ 

j-g.  =  cos  (Pc  -  (ut  -  Uc)  jy- 


(21) 


(22) 


Let  the  projection  of  rB  (figure  2(a))  on  a  radial  line  be  u^'(^).  Then,  provided 
there  is  neither  clearance  nor  interference  between  the  tube  and  die  entrance  guide 
and  the  dimensionless  die  radius  R  is  constant,  the  die  angle  i  must  satisfy 


or 


u^'  «  R(1  -  cos  >]/’) 


>|(  =  2  arcsin 


(23) 


From  geometry  u^*  depends  on  u^  and  cp^, 
*  u^  -f  h(l-cos 


(2U) 
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The  ratio  of  shearing  force  to  bending  moment  at  the  crack  tip  depends  on  the  geometry 
of  the  petal  and  die  and  the  friction  angle  (figure  2(b)).  This  ratio  becomes 


(25) 

(26) 

Finally,  one  obtains  another  expression  relating  b  and  p  from  equations  (80)  and 
(8l)  of  the  Appendix  with  r}  =  p. 

(27) 

"2  l-p2-|l/3p3 

Fourteen  variables  U,  p,  q,  f,  u^,  u^,  f,  v,  and  b)  and 

thirteen  functional  relationships  [equations  (l),  (12)  through  (15),  and  (20) 
through  (2?)]  have  been  introduced.  The  constant  parameters  of  the  problem  are 
R,  h,  u^,  and  [i.  The  functions  <Pq(^)  a.nd  specify  the  initial  state.  The 

seven  variables  p,  x^,  f,  v,  u^*,  and  b  may  be  eliminated  by  algebraic  processes, 
leaving  six  coupled  functional  relationships  among  the  remaining  variables. 

Solution 


2b  = 


m. 


cos  (v  -  (pc) 

Xz^Hr"h  sin  <p^)cos  v  +  (u^  -  u^^  -  h  cos  (p^)  sin  v 


where 


V  =  \(r  +  tan”^  M- 


It  is  most  convenient  to  treat  ^  as  the  independent  variable.  A  numerical 
step-by-step  method  of  solution,  beginning  at  0  and  advancing  by  small  steps 

has  been  developed.  The  method  employs  the  trapezoid  rule  of  integration  and 
reiteration  of  each  step  until  successive  values  of  all  variables  agree  within  a 
specified  percentage  error.  It  is  desirable  to  decrease  the  size  of  the  step  LI 
as  the  integration  proceeds  because  as  p  becomes  small  f  becomes  large.  In  the 
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numerical  procedure  an  initial  value  of  is  chosen,  and  LI  is  halved  whenever 

the  value  of  f  increases  by  more  than  a  specified  percentage  over  one  step.  As  p 
approaches  zero,  the  behavior  becomes  singular;  therefore,  the  procedure  must  be 
terminated  when  p  becomes  small.  Then  an  asymptotic  solution  employing  a  variable 
other  than  t  as  the  independent  variable  is  used.  This  asymptotic  method  is 
described  later. 

The  method  was  programed  in  the  Fortran  language  for  an  IBM  7070  computer. 

After  a  few  trials,  it  was  found  that  satisfactory  results  were  obtained  by  choosing 
LI  =  0.1  initially,  halving  LI  when  f  changed  by  more  than  2$%y  and  reiterating  to 
an  accuracy  of  0.3^.  The  procedure  was  terminated  when  either  p  became  less  than 
,03  or  became  less  than  ,0005*  The  functions  chosen  for  (p^  and  Xq  are  the  dimen¬ 
sionless  slope  and  curvature  resulting  from  the  application  of  a  transverse  shear 
force  at  the  end  of  the  tube  to  open  it  up  by  a  radial  amount  u^. 


Asymptotic  Method  for  Small  p 

From  the  differentiation  of  (13)  and  (Ih)  with  respect  to  ^  and  the  condition 
that  p  is  small  one  finds 

K  d'f’o  ‘Pc 

TT  *  p 

'^’‘c  ~  2  rr'Pc 

~  TT  3"^^  P 

In  general,  where  the  asymptotic  solution  applies,  and  vanish.  Let  this 

be  the  case.  Now  with  as  the  independent  variable, 


dx 


c 


d(p 


c 


2 

3 


(28) 


(29) 


ii.  ~  E_ 

<iPc  ~  <Pc 

By  changing  the  indepencient  variable  from  ^  to  cp^  in  (20),  (21),  and  (22)  and 
employing  (29)  one  finds 


dP 

c 


:  1 


(30) 


du. 

X»  «r 


(31) 


~  -  Up 

c 


For  small  p  equation  (27)  reduces  to 


(32) 


•i  ^ 

b  A*  p 


(33) 


One  must  now  solve  the  eight  coupled  functional  relationships  (23)  through  (26) 
and  (30)  through  (33)  for  eight  dependent  variables  (p^,  u^,  z^,  p,  b,  u^', 
and  V  as  functions  of  the  independent  variable  <P^,  Again  a  numerical  step-by- 
step  procedure  ett^loying  the  trapezoid  integration  rule  and  reiteration  is  used* 

The  initial  data  for  this  computation  is,  of  course,  taken  from  the  final  step  of 
the  previous  computation.  A  step  of  A<p^  »  0.05  was  used.  The  procedure  is  repeated 
until  p  becomes  zero.  Both  ^  and  may  be  calculated  afterward  from  (28)  and 
(29).  Equation  (29)  shows  that  the  cracks  are  arrested  when  p  reaches  zero,  i.e., 
the  crack  tips  and  plastic  hinges  have  the  same  C -coordinate  and  the  shearing  force 
at  the  petal  bases  vanishes. 
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Ill  Petal  Bending  with  Stationary  Cracks 

After  the  crack  has  been  arrested,  the  plastic  hinge  moves  away  from  the  crack 
tip  into  the  petal.  Figure  U(a)  shows  a  cross  section  of  the  petal  with  the  plastic 
hinge  at  point  H,  the  crack  tip  being  at  point  C.  It  is  assumed  that  the  petal 
behaves  as  a  plastic  curved  beam  subjected  to  tip  loading.  The  bending  moment  is 
a  maximum  at  €  =  g,  the  location  of  the  plastic  hinge.  The  shear  force  vanishes 
at  ?  =  g;  hence,  from  figure  U(b) 

>1^  +  tan”^  M--  <P.  ■  J  .  (3U) 


For  a  constant  friction  coefficient 


d 


A  change  in  the  slope  angle  at  the  plastic  hinge  is  the  sum  of  a  rigid 

body  rotation  Acp.  of  the  region  between  the  hinge  and  the  petal  tip  and  a  change 
t 

-xAg  due  to  the  movement  Ag  of  plastic  hinge  along  the  petal  with  curvature  x  . 
In  derivative  form  this  condition  becomes 


Other  geometrical  relationships  are 


^  d<p^ 

d<^  -y 


^._||_C0S 


(39) 


sin  (p^ 


where  y  is  the  axial  projection  of  line  TH  of  figure  U(a)  and  is  the  total 

radial  displacement  at  the  plastic  hinge.  The  die  contact  angle  \|^  is  related  to 
the  radial  displacement  at  the  petal  tip  by  equations  (23)  and  (2li).  Now  substituting 
(2U)  into  (23),  differentiating  with  respect  to  (p^,  using  the  result  of  (35),  and 
rearranging,  one  finds 


R  sin  i 
y  +  h  sin 


(ho) 


Then  from  (36)  and  (UO) 


i_  [1  _  R  sin  t  , 

X  y  +  h  sin  <p^ 


(lii) 


In  (Ul)  X  is  the  curvature  of  the  petal  just  ahead  of  the  moving  hinge.  Its  value 
is  obtained  from  the  solution  of  the  crack  growth  problem. 

The  system  of  six  differential  equations  (35)  and  (37)  through  (Ul)  can  be 
solved  for  the  six  variables  t,  g,  y,  and  u^  as  functions  of  (p^.  The 
trapezoid  rule  of  numerical  integration  and  reiteration  of  each  step  to  a  desired 
accuracy  have  been  used  to  obtain  step-by-step  numerical  solutions.  The  initial 
data  used  for  these  calculations  is  the  terminal  data  taken  from  the  s  olution  of  the 


crack  growth  problem, 


IV  Load-Displacement  Behavior 


Since  it  has  been  assumed  that  contact  between  the  tube  material  and  the  die 
occurs  only  at  the  petal  tips,  the  axial  force  per  unit  circumference  in  the  rigid 
tube  depends  only  on  the  magnitude  of  the  die  contact  force  F  per  unit  circum¬ 
ference  of  undeformed  tube  and  the  angle  v*  The  compressive  stress  in  the  tube  is 


sin  V  • 


(1*2) 


For  the  crack  propogation  problem  one  finds  F  by  eliminating  from  the 

expressions 


4m„  =  1  -  p2  -  ^/3  p^ 


(U3) 


and 


^  [(z^  +  h  sin  tp^)  cos  V  +  (u^  -  -  h  cos  cp^)  sin  v]  (Ul*) 


a  t2 


or  by  eliminating  q  ’from  the  expressions 
c 


2q^-p  *  -y  VJ  p 


4-  ^/3  2 


(U5) 


and 


_  — cos  (v  -  «p_). 


a 


(1*6) 


After  the  crack  has  been  arrested,  the  condition  that  ^ini^  ■  1  yields  the  result 


F  n/aT 


0.20 


(1*7) 


2  (y  +  h  sin  <p.)  cos  v  +  (u  -u  -  h  cos  cp  )  sin  v 

t/  li  tl  w 
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The  meridianal  membrane  stress  at  the  petal  root  is 


o 

The  yield  condition  used  so  far  is  independent  of  this  stress.  Such  a  condition  is 
valid  provided  |n  |  is  small  in  comparison  with  one.  For  the  present  problem 
this  is  not  always  the  case;  therefore,  a  correction  for  this  effect  is  needed.  To 
obtain  an  approximate  correction,  assume  that  the  effect  of  n^  on  the  projection 
of  the  yield  surface  on  the  m,  n^  plane  is  to  change  only  the  size  but  not  the 
shape  of  the  figure.  The  plastic  analysis  of  such  a  cylindrical  shell  is  given  in 
the  Appendix.  The  intersection  of  the  m-axis  and  the  yield  surface  is  given  by 

2 

Urn  =  1  -  n 

X 


p 

hence,  the  size  factor  for  the  yield  surface  is  assumed  to  be  1  -  n^  •  For  such 

a  yield  surface  m  and  q  of  equations  (69)  and  (70)  are  replaced  by  m/(l  -  n  ) 
2 

and  q/(l  -  n^  ),  respectively.  To  employ  such  a  yield  surface  in  the  crack 

p 

propagation  analysis,  replace  m  in  equation  (h3)  by  m  /(l  -  n  and  q  in 

c  c  xc  c 

2 

equation  (U5)  by  q  /(l  -  n  ),  The  above  approximate  yield  surface  is  expected 

to  give  results  very  close  to  those  obtained  from  the  exact  Tresca  yield  surface 

because  the  effect  of  I n  I  on  the  yield  surface  becomes  appreciable  only  when 

xc 

the  generalized  stress  profile  approaches  the  m-axis,  where  the  exact  and  approx¬ 
imate  yield  surfaces  coincide.  To  apply  the  approximate  correction  for  axial 
force  to  the  stationary  crack  problem,  multiply  the  numerator  on  the  right  hand  side 
of  equation  (U7)  by  (l  -  )  where 


(li9) 
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The  coti5)ressive  stress  in  the  tube  found  by  employing  the  axial  force  dependent 
yield  surface  is  denoted  by  o'. 

Let  D  denote  the  dimensionless  displacement  of  the  undeformed  part  of  the 
tube  relative  to  the  die  and  take  D  =  0  when  the  tip  enters  the  curved  part  of 
the  die.  It  is  assumed  that  there  is  neither  clearance  nor  interference  between 
the  tube  and  the  die, entrance  and  that  meridianal  lines  on  the  tube  middle  surface 
do  not  change  in  length.  From  geometry  of  the  crack  growth  problem 

D  =  ^  _  z  -  h  sin  <p.  +  R  sin  \|f  . 
b  w 

Forttie  problem  of  petal  bending  with  stationary  cracks  the  derivative  of  D  with 
respect  to  the  independent  variable  is 


dP 

dPh 


=  R  cos  +  (u^ 


d(Pt 

%  -  ^  ‘Pt^  ^ 


(51) 


The  rigid  tube  displacement  relative  to  the  die  is  then  found  by  numerical 
integration. 


17 


V  A  Theoretical  Analysis  of  the  Rolling  Process 


For  small  t/fe  ratios  the  a  priori  condition  h  does  not  hold.  In  this 

case  the  petals  formed  are  likely  to  roll  up  to  a  radius  close  to  the  die  radius. 

This  phenomenon  is  called  rolling.  In  the  following  analysis  it  is  assumed  that, 

when  rolling  occurs,  the  flow  of  material  in  the  deforming  region  is  steady  and 

that  the  curvature  of  the  formed  petal  corresponds  to  the  die  radius.  The 

satisfaction  of  both  of  these  conditions  places  some  restrictions  upon  the  die 

configuration  near  the  entrance.  This  restriction  will  be  discussed  later.  Figure  5 

shows  the  petals  bent  to  the  die  radius.  The  plastic  analysis  of  cylindrical  shells 

of  the  Appendix  applies  to  the  deforming  portion  of  the  tube.  The  tub'e  material 

flows  steadily  along  a  path,  entering  the  deforming  region  at  t)  =*  0  and  leaving 

at  Tl  “  'H  ^  the  coordinate  at  the  crack  tip  or  petal  base.  Here  "H  is  the 
c 

dimensionless  coordiriate  of  a  point  of  the  deforming  tube  measured  in  a  frame  of 
reference  fixed  to  the  die.  The  point  at  which  plastic  deformation  commences  is  the 
origin.  When  axial  compression  is  neglected,  the  average  axial  velocity  over  the 
cross  section  is  independent  of  t]  ,  The  time  required  for  a  material  point 
entering  the  deforming  region  at  t]  =  o  to  reach  the  position  'H  is  proportional 
to  T]  .  Then 

TT  =  (52) 

w 

where  w  is  the  dimensionless  velocity  of  the  tube  relative  to  the  die.  The 
dimensionless  radial  displacement  u  of  a  point  at  r)  is  found  by  integrating 
the  radial  velocity  v  with  respect  to  time.  The  dependence  of  the  radial 
velocity  on  the  coordinate  t]  is  given  by  equation  (77)  or  (83)  in  the  Appendix. 

For  steady  flow  A  must  be  independent  of  .  The  expression  for  the  radial 
displacement  becomes 
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u  (ti)  =  X'  (4  Ti^  +  ^  '^)  H  (ti)  i^h) 

where  \*  «  A/w.  Expressions  for  the  slope  and  the  curvature  are  found  by 
differentiation  with  respect  to  t)  .  Thus 

(p  (p)  «  X*  (t)  +  4  H  (t))  (55) 

X  (t))  =  X‘  (l  +  t)  +  7]^)  H  (p)  .  (56) 


Two  conditions  are  maintained  at  P  =  'P  .  They  are  that  the  petal  curvature 

c 

corresponds  to  the  die  radius  and  that  the  radial  displacement  is  that  which  is 
required  to  maintain  the  crack.  These  two  conditions  are  expressed  by 


(1  -K 


1 

“T 


(57) 


1  : 


)  =  ac 


n  and  X'  are  found  by  solving  equations  (5?)  and  (58)  simultaneously. 

The  slope  of  both  the  die  and  the  petal  must  be  the  same  at  point  C  if  the 
petal  is  to  match  the  die  exactly.  This  cannot  be  the  case  if  the  die  radius  is 
constant.  If  the  die  radius  is  modified  along  the  arc  BC  as  shown  in  Figure  6, 
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the  condition  of  matching  slopes  is  satisfied.  _For  slopes  satisfying  cp^  «  1 


“o  ■  i 


from  which 

R  -  h  2  u  X 
1  c  c 

-  . .  S  I  |» 

R  -  h  <P  ^ 

c 


($9) 


The  numerical  value  of  (R^  -  h)/(R  -  h)  lies  between  1.0  and  1.2j  therefore, 
dies  by  which  steady  state  formation  of  petals  of  constant  curvature  may  be  achieved 
are  almost  constant  radius  dies. 

Figure  7  shows  a  free  body  diagram  of  the  petal  which  may  carry  an  arbitrary 
distribution  of  contact  force  p  and  a  friction  force  [xp.  If  there  is  no  friction, 
the  tangential  force  must  vanish.  For  the  frictionless  die 


N 

c 


Rp  sin  (<P  -  cp  )  d(p 
c 


In  terms  of  dimensionless  variables 


n 

c 


m 


R  -  h 


(60) 


where  m  is  found  from  equation  (69)  of  the  Appendix,  with  tj  =  t)  .  The  stress 
c  c 

in  the  undeformed  tube  is 


n  cos  <p  +  — r—  sin  (p 

C  C  U  c 


(61) 


where  is  found  from  equation  (70)  of  the  Appendix. 


With  some  algebraic 
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manipulation  equation  (6l)  reduces  to 


1 

TT 


■n  6  . 

C  C 


(62) 


For  practical  values  of  t)^  the  factor  in  parentheses  in  (62)  has  a  value  in  the 

range  1.0  to  1.33.  A  limitation  of  the  foregoing  analysis  is  that  must 

satisfy  e  4  t/(R  -  h). 
c 

For  6  >  i  t/(R  -  h)a  portion  of  the  unsplit  but  deformed  tube  has  contact  with 

c 

the  die.  Again  shell  theory  is  used  to  determine  the  stress  in  the  tube.  The 
resulting  expression  for  the  stress  in  the  tube  is 


^  .  1 
o 


(63) 


The  details  of  this  calculation  are  omitted  here  for  brevity. 

When  a  small  coefficient  of  friction  (about  *2  or  less)  between  the  die  and 
petal  exists,  equation  (60)  must  be  replaced  by 


Equation(6U)  is  not  an  exact  expression,but  contains  only  the  first  order  terms  in  |ji„ 
The  derivation  of  equation  (6U)  is  omitted  for  brevity.  Figure  8  shows  the  results 
of  calculations* 
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71  Theoretical  Results  and  Comparison  with  Experiment 


Geometric  results  of  a  sample  calculation  are  shown  in  Figure  9.  The  data 

a/R  =  1.7$,  t/R  =  .6,  and  e  =  .02  are  typical  values  for  the  tube  fragments- 

c 

tion  tests  of  reference  [2]  using  202U  -  T3  aluminum  alloy  as  the  tubing 

material.  The  value  for  e  was  obtained  from  measurements  of  deformed  specimens 

c 

of  these  tests.  The  friction  coefficient  is  assumed  to  vanish  corresponding  to  a 
well  lubricated  die.  Figure  10  shows  how  the  C -coordinate  of  the  plastic  hinge  g 
and  that  of  the  crack  tip  ^  change  as  the  tube  is  pushed  over  the  die.  Figure  10 
also  shows  the  load- displacement  behavior  of  the  rigid  tube.  In  Figures  9  and  10  the 
dotted  lines  at  the  left  show  interpolations  to  the  known  initial  data  for  D  =  0, 

The  interpolation  line  for  vs,  D  in  Figure  9  begins  at  the  value  of  D  for 

which  g  reaches  a  maximum  as  would  be  expected.  The  area  under  the  <7/7^  vs,  D 
curve  in  Figure  10  represents  the  energy  absorbed  during  the  fragmentation  process. 
The  maximum  ordinate  of  the  vs,  D  curve  indicates  the  maximum  compressive 

stress  reached  in  the  undeformed  part  of  the  tube.  Buckling  of  the  tube  as  a  short 
column  may  occur  if  becomes  too  large.  For  the  buckling  problem  it  is 

necessary  to  treat  the  tube  material  as  an  elastic-work  hardening  material.  The 
buckling  problem  is  not  considered  in  this  paper. 

Figures  11(a)  and  11(b)  show  the  petal  shapes  at  the  time  when  the  crack  becomes 
arrested  and  when  the  motion  of  the  plastic  Viinge  in  the  petal  changes  its  direction. 
These  figures  show  the  highly  non-uniform  distribution  of  bending  in  the  petals 
for  this  case.  Transverse  fracture  of  the  petal  is  very  likely  when  the  configura¬ 
tion  of  Figure  11(b)  is  reached  due  to  the  large  bending  strains  developed. 

Figure  12  shows  the  effect  of  the  ratio  t/R  on  the  load-displacement 
behavior  of  tubes  pushed  over  a  die.  For  the  three  cases  shown,  the  constant  data 

is  Sl/R  =  1,75j  e  =  ,02,  and  p,  =  0,  For  the  cases  t/R  =  ,3  and  t/R  =  ,U6, 
c 

condition  h  (that  contact  between  the  tube  material  and  the  die  occurs  only  at 
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the  tips  of  the  fortned  petals)  is  valid  only  during  the  early  stages  of  petal  growth. 
The  load-displacement  behavior  for  these  cases  is  shown  only  up  to  the  displacement  for 
which  condition  h  holds.  At  larger  displacements  contact  occurs  at  both  the  petal 
tip  and  over  a  circle  in  the  unsplit  portion  of  the  tube.  Continued  displacement  of 
the  tube  would,  due  to  the  new  contact  force,  cause  an  increase  of  the  shearing 
force  on  the  unsplit  portion  of  the  shell.  This  in  turn  causes  an  increase  of 
length  of  the  plastic  region  in  the  unsplit  portion  of  the  tube.  It  is  conjectured 
that  then  the  length  of  the  plastic  region  in  the  unsplit  portion  of  the  tube  would 
reach  a  relative  maximum  and  then  decrease  due  to  rapid  running  of  the  meridianal 
cracks.  A  theoretical  analysis  of  this  deformation  would  be  too  complicated  to  carry 
out5  however,  reasonable  estimates  of  the  force  level  in  the  undeformed  tube  may  be 
obtained  by  supposing  that  the  length  of  plastic  region  remains  constant,  i.e.,  by 
using  the  results  for  steady  state  rolling. 

Figure  13  shows  how  a  change  in  the  ratio  a/R  affects  the  load-displacement 

behavior  for  t/R  =  .U6,  =  .02,  and  p.  =  0.  Figures  lit  and  IS  show  how  changing 

the  friction  coefficient  from  0  to  0,2  affects  the  load-displacement  behavior 

when  6  =  .02,  a^  =  1,75,  and  t/R  =  ,U6  or  t/R  =  .6.  Note  the  large  effect  of 

c 

friction  when  t/R  =  ,6. 

Results  of  tests  of  tubes  of  202U  -T3  aluminum  alloy  have  been  communicated 
to  the  author  [2 ].  Figure  16  shows  the  stress-strain  behavior  in  conpression  for 
this  material  and  that  for  a  rigid  perfectly  plastic  idealization  with  a  flow  stress 
of  60  ksi.  The  dotted  lines  of  Figure  17  show  the  measured  stress  levels  plotted 
as  a  function  of  t/(R  -  h)  for  various  values  of  a/(R  -  h).  These  curves  were 
plotted  from  an  empirical  expression  given  in  [2]  by  taking  the  flow  stress 
as  60,000  psi.  The  slight  increase  of  with  increasing  a/(R  -  h)  is  in 

agreement  with  the  theoretical  results  of  Figure  13.  Curve  1  of  Figure  17  shows 
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theoretical  stress  levels  for  steady  state  rolling  and  curve  2  shows  that  for 
fragmentation.  For  ,35  <  t/(R  -  h)  <  ,7  agreement  between  the  theoretical 
fragmentation  stress  and  the  measured  stress  is  goodo  For  larger  t/(R  -  h)  the 
theory  predicts  a  stress  in  the  undeformed  part  of  the  tube  which  is  much  less  than 
the  measured  values. 

Some  reasons  for  the  difference  are  given  in  the  following: 

1,  For  large  t/(R  -  h)  the  theoretical  analysis  shows  that  the  contact 
force  at  the  petal  tips  becomes  large,  so  that  the  edges  would  be  crushed.  The 
energy  absorbed  by  local  plastic  flow  and  the  change  in  petal  shape  due  to  this 
crushing  have  not  been  considered  in  the  theoretical  analysis. 

2,  The  cross  section  of  the  petals  is  not  rectangular  in  shape  but  is  a  sector 
of  an  annulus  which  has  a  greater  plastic  limit  moment  than  the  assumed  rectangular 
section, 

3,  The  dies  were  machined  to  have  a  slight  interference  with  the  inside  of  the 
tubular  specimens  along  the  die  entrance  guide.  Considerable  friction  may  occur  along 
this  contact  surface. 

lio  The  assumption  that  the  length  of  the  middle  surface  does  not  change  is 
poorest  for  the  case  of  large  t/R;  hence,  the  calculated  values  of  D  are  too 
small. 
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VII  Conclusions 


1,  In  the  fragmentation  regiite  meridianally  running  cracks  advance  rapidly  when  the 
shearing  force  at  the  petal  roots  is  large  and  become  arrested  when  this  shear- 
ing  force  vanishes, 

2,  Large  localized  bending  strains  develop;  therefore,  points  where  transverse 
fracture  is  likely  to  occur  may  be  determined, 

3,  The  direct  stress  in  the  tube  increases  with  increasing  t/R  in  both  the  frag¬ 
mentation  and  the  rolling  regimes, 

Uo  The  direct  stress  in  the  tube  increases  slightly  with  increasing  a/R, 

5,  The  effect  of  friction  on  the  force  level  in  the  undeformed  part  of  the  tube  is 
significant, 

6,  The  compressive  stress  in  the  undeformed  part  of  the  tube  may  reach  peak  values 
sufficiently  large  to  cause  column  buckling  of  the  tube  unless  the  tube  is 
supported  laterally  to  prevent  such  buckling, 

7,  For  large  t/R  the  energy  absorbed  by  crushing  at  the  petal  tips  may  be  signi¬ 
ficant. 
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APPENDIX  Plastic  Analysis  of  a  Cylindrical  Shell 


A  semi -infinite  cylindrical  shell  is  subjected  to  the  axisymmetric  edge  bending 
moment  M  and  shearing  force  Q  as  shown  in  Figure  18(a),  The  coordinate  x  is 
measured  from  the  point  where  the  shearing  force  vanishes.  The  following  dimension¬ 
less  notation  is  introduced: 


m  =  M/gt  T^  , 
o  ^ 


q  “  Q  n/at/o  T' 


n9=NeAT  , 


n  =  N  /o-  T 

X  X  o 


•n  =  x/n/at 


In  terms  of  the  dimensionless  quantities ^  the  equations  of  equilibrium  are 


When  n^  =  0^  ni  >  0,  nQ>  the  Tresca  yield  condition  [Figure  18(f)]  for  a 


cylindrical  shell  is  (see  reference  [i3]) 


^  (m,  xiq)  =  0 


where 


$  =  Urn  +  (2n  -  1)"^  -  1, 
9 
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The  solution  to  (65),  (66),  and  (6?)  is 


Urn  =  1  -  T)^ 


(69) 


2q  =  T)  +  I  r)^  +  T}^  (70) 

^0  =  I  +  In  +  •  (71) 


Figures  18(b),  l8(c),  and  l8(d)  show  the  variations  of  m,  q,  and  with  t].  By 
eliminating  t)  from  equation  (69)  and  (70)  one  obtains  the  first  quadrant  of  the 
yield  curve  in  load  space  shown  in  Figure  18(e). 

For  a  cylindrical  shell  the  generalized  strain  rates  are  the  curvature  rate  K 

* 

and  the  circumferential  strain  rate  e  •  These  generalized  strain  rates  are  related 
to  the  radial  velocity  V  by 


d^  V 

dx^ 


1 

AT 


*1^ 


(72) 


e 


(73) 


The  flow  rule  of  the  theory  of  plasticity  states  that  the  generalized  strain  rate 
vector  points  in  the  direction  of  the  outward  normal  to  the  yield  surface  at  a 
regular  point.  This  means  that 


a 

o 

1  ^  X  0^ 

where  X  >  0* 


(7U) 

(75) 
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Expression  (68)  is  substituted  into  (7U)  snd  (7^)  snd  is  eliminated  to 
obtain  a  differential  equation  for  V. 


(76) 


Let  V  =  V/ be  the  dimensionless  radial  velocity,  w  be  the  dimensionless 

•  -■  ■■  • 

angular  velocity  and  "x-  =vAT  K  be  the  dimensionless  curvature  rate.  The  solution 
to  (76)  is 


V  =  A(ti  +  T)^)  H  (r)) 


(77) 


03  =  =  A(1  +T1  +  -5-11^)  H  (ti) 


(78) 


;  =  -|^  =  a[6(ti)  +  (1  +  ^-n)  H  (n)]  (79) 

where  A  >  0^  H  (t}  )  is  the  Heaviside  step  function  and  6  (t]  )  is  the  Dirac  delta 
function. 

Equation  (71)  shows  that  the  variation  of  3:^  with  t]  is  almost  linear.  To 
obtain  suitable  approximations  by  polynomials  of  a  lower  order,  assume  that  the 
variation  of  nQ  with  t]  is  linear,  equations  (65)  and  (66)  hold,  and  that 
^0  =  4  where  m  ==  J  and  =  1  where  m  =  0.  The  resulting  approximations  are 


(80) 

n  +  ^fT r\^ 

(81) 

1  +  T]  . 

(82) 
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An  approximate  expression  for  the  velocity  field  is  found  by  observing  the 
similarity  of  the  exact  expressions  for  v  and  2q  and  applying  this  similarity  to 
the  approximate  expression.  Thus  the  approximate  velocity  field  is 


=  A  (t)  +  H  (ri) 

(83) 

-  A  (1  1  Ni3  Tl)  H  (t)) 

(8U) 

=  A  [6(ti)  +  |Nf3  H  (Tl)] 

(85) 

When  n^  is  not  zero  the  yield  surface  of  a  cylindrical  shell  of  Tresca  material 
subjected  to  axisymmetric  loading  has  the  form 

(m,  Hq,  n^)  =  0 

An  approximate  expression  for  p  suitable  for  the  region  Uq  >  0,  m  >  0, 

n  X  <  0  small,  is 


is  obtained  by  replacing  m  and  nn  in  the  expression  ^  by  - -  and  - r*  , 

1  1  -  V 

respectively.  The  projection  of  =  0  on  the  plane  n^  =  -0,5  is  shown  in  figure  18(f) , 

The  plastic  analysis  involves  the  satifaction  of  equations  (65), (66),  (6?),  and  (86), 

The  solution  is  given  by  expressions  (69),  (70),  and  (71)  with  m,  q,  and  replaced 

riQ 

■by -  ^  - 3 -  ^  and -  ,  respectively.  Fig.  I8  (e)  shows  such  a  yield 

1  -  n^  1  -  1  - 

curve  in  load  space  for  =  -.5.  For  the  approximate  yield  surface,  the  normality 
condition  yields  the  transverse  velocity  field  given  by  equations  (77),  (78),  and 
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(79)  with  no  changes. 

The  derivations  and  results  of  this  Appendix  were  guided  by  [b]  . 
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CRACK  TIP 


PLASTIC  HINGE 


PLASTIC 

REGION 


PLASTIC  HINGE 


CRACK  TIP 


Fig.  5 


Cross  section  of  tube,  petal,  and  die  during  rolling 


Fig.  11  Computed  petal  shapes:  t/R  =  .6,  a/R  -  1.75>  ^.nd  p.  -  0. 

(a)  'When  the  cracks  are  arrested^  (b)  When  the  plastic  hinges  stop  moving 


Effect  of  friction  on  tube  stress-tube  displacement  behavior:  t/R 
a/R  =  1.75,  and  e  =  ,02. 


b  b°  b  I  b' 


Fig*  l8  (a)  Shell  and  loading^  (b)  dimensionless  bending  moment^  (c)  dimen 
sionless  shearing  force^  (d)  dimensionless  circumferential  tension^ 
(e)  yield  condition  in  load  space,  (f)  Tresca  yield  condition  for  a 
cylindrical  s hello 


NASA- Langley,  1966  CR-442 
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